Molecular mechanism of dynein-dynactin complex assembly by LIS1.

Cytoplasmic dynein is a microtubule motor vital for cellular organization and division. It functions as a ~4-megadalton complex containing its cofactor dynactin and a cargo-specific coiled-coil adaptor. However, how dynein and dynactin recognize diverse adaptors, how they interact with each other during complex formation, and the role of critical regulators such as lissencephaly-1 (LIS1) protein (LIS1) remain unclear. In this study, we determined the cryo-electron microscopy structure of dynein-dynactin on microtubules with LIS1 and the lysosomal adaptor JIP3. This structure reveals the molecular basis of interactions occurring during dynein activation. We show how JIP3 activates dynein despite its atypical architecture. Unexpectedly, LIS1 binds dynactin's p150 subunit, tethering it along the length of dynein. Our data suggest that LIS1 and p150 constrain dynein-dynactin to ensure efficient complex formation.

Rapid Optogenetic Clustering in the Cytoplasm with BcLOVclust.

Protein clustering is a powerful form of optogenetic control, yet remarkably few proteins are known to oligomerize with light. Recently, the photoreceptor BcLOV4 was found to form protein clusters in mammalian cells in response to blue light, although clustering coincided with its translocation to the plasma membrane, potentially constraining its application as an optogenetic clustering module. Herein we identify key amino acids that couple BcLOV4 clustering to membrane binding, allowing us to engineer a variant that clusters in the cytoplasm and does not associate with the membrane in response to blue light. This variant-called BcLOVclust-clustered over many cycles with substantially faster clustering and de-clustering kinetics compared to the widely used optogenetic clustering protein Cry2. The magnitude of clustering could be strengthened by appending an intrinsically disordered region from the fused in sarcoma (FUS) protein, or by selecting the appropriate fluorescent protein to which it was fused. Like wt BcLOV4, BcLOVclust activity was sensitive to temperature: light-induced clusters spontaneously dissolved at a rate that increased with temperature despite constant illumination. At low temperatures, BcLOVclust and Cry2 could be multiplexed in the same cells, allowing light control of independent protein condensates. BcLOVclust could also be applied to control signaling proteins and stress granules in mammalian cells. While its usage is currently best suited in cells and organisms that can be cultured below ∼30 °C, a deeper understanding of BcLOVclust thermal response will further enable its use at physiological mammalian temperatures.

Autoinhibition in the Signal Transducer CIN85 Modulates B Cell Activation.

Signal transduction by the ligated B cell antigen receptor (BCR) depends on the preorganization of its intracellular components, such as the effector proteins SLP65 and CIN85 within phase-separated condensates. These liquid-like condensates are based on the interaction between three Src homology 3 (SH3) domains and the corresponding proline-rich recognition motifs (PRM) in CIN85 and SLP65, respectively. However, detailed information on the protein conformation and how it impacts the capability of SLP65/CIN85 condensates to orchestrate BCR signal transduction is still lacking. This study identifies a hitherto unknown intramolecular SH3:PRM interaction between the C-terminal SH3 domain (SH3C) of CIN85 and an adjacent PRM. We used high-resolution nuclear magnetic resonance (NMR) experiments to study the flexible linker region containing the PRM and determined the extent of the interaction in multidomain constructs of the protein. Moreover, we observed that the phosphorylation of a serine residue located in the immediate vicinity of the PRM regulates this intramolecular interaction. This allows for a dynamic modulation of CIN85's valency toward SLP65. B cell culture experiments further revealed that the PRM/SH3C interaction is crucial for maintaining the physiological level of SLP65/CIN85 condensate formation, activation-induced membrane recruitment of CIN85, and subsequent mobilization of Ca2+. Our findings therefore suggest that the intramolecular interaction with the adjacent disordered linker is effective in modulating CIN85's valency both in vitro and in vivo. This therefore constitutes a powerful way for the modulation of SLP65/CIN85 condensate formation and subsequent B cell signaling processes within the cell.

Integrated virtual screening and MD simulation approaches toward discovering potential inhibitors for targeting BRPF1 bromodomain in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most aggressive and life-threatening cancers. Although multiple treatment options are available, the prognosis of HCC patients is poor due to metastasis and drug resistance. Hence, discovering novel targets is essential for better therapeutic development for HCC. In this study, we used the cancer genome atlas (TCGA) dataset to analyze the expression of bromodomain-containing proteins in HCC, as bromodomains are emerging attractive therapeutic targets. Our analysis identified BRPF1 as the most highly upregulated gene in HCC among the 43 bromodomain-containing genes. Upregulation of BRPF1 was significantly associated with poorer patient survival. Therefore, targeting BRPF1 may be an approach for HCC treatment. Previously, several potential inhibitors of BRPF1 bromodomain have been discovered. However, due to the limited clinical success of the current inhibitors, we aim to search for new inhibitors with high affinity and specificity for the BRPF1 bromodomain. In this study, we utilized high-throughput virtual screening methods to screen synthetic and natural compound databases against the BRPF1 bromodomain. In addition, we used machine learning-based QSAR modeling to predict the IC50 values of the selected BRPF1 bromodomain inhibitors. Extensive MD simulations were used to calculate the binding free energies of BRPF1 bromodomain and inhibitor complexes. Using this approach, we identified four lead scaffolds with a similar or better binding affinity towards the BRPF1 bromodomain than the previously reported inhibitors. Overall, this study discovered some promising compounds that have the potential to act as potent BRPF1 bromodomain inhibitors.

Mechanistic Insights into the Interactions of Arl8b with the RUN Domains of PLEKHM1 and SKIP.

Arl8b, a specific Arf-like family GTPase present on lysosome, and plays critical roles in many lysosome-related cellular processes such as autophagy. The active Arl8b can be specifically recognized by the RUN domains of two Arl8b-effectors PLEKHM1 and SKIP, thereby regulating the autophagosome/lysosome membrane fusion and the intracellular lysosome positioning, respectively. However, the mechanistic bases underlying the interactions of Arl8b with the RUN domains of PLEKHM1 and SKIP remain elusive. Here, we report the two high-resolution crystal structures of the active Arl8b in complex with the RUN domains of PLEKHM1 and SKIP. In addition to elucidating the detailed molecular mechanism governing the specific interactions of the active Arl8b with the RUN domains of PLEKHM1 and SKIP, the determined complex structures also reveal a general binding mode shared by the PLEKHM1 and SKIP RUN domains for interacting with the active Arl8b. Furthermore, we uncovered a competitive relationship between the RUN domains of PLEKHM1 and SKIP in binding to the active Arl8b as well as a unique small GTPase-binding mode adopted by the PLEKHM1 and SKIP RUN domains, thereby enriching the repertoire of the RUN domain/small GTPase interaction modes. In all, our findings provide new mechanistic insights into the interactions of the active Arl8b with PLEKHM1 and SKIP, and are valuable for further understanding the working modes of these proteins in relevant cellular processes.

Structural insights of the p97/VCP AAA+ ATPase: How adapter interactions coordinate diverse cellular functionality.

p97/valosin-containing protein is an essential eukaryotic AAA+ ATPase with diverse functions including protein homeostasis, membrane remodeling, and chromatin regulation. Dysregulation of p97 function causes severe neurodegenerative disease and is associated with cancer, making this protein a significant therapeutic target. p97 extracts polypeptide substrates from macromolecular assemblies by hydrolysis-driven translocation through its central pore. Growing evidence indicates that this activity is highly coordinated by "adapter" partner proteins, of which more than 30 have been identified and are commonly described to facilitate translocation through substrate recruitment or modification. In so doing, these adapters enable critical p97-dependent functions such as extraction of misfolded proteins from the endoplasmic reticulum or mitochondria, and are likely the reason for the extreme functional diversity of p97 relative to other AAA+ translocases. Here, we review the known functions of adapter proteins and highlight recent structural and biochemical advances that have begun to reveal the diverse molecular bases for adapter-mediated regulation of p97 function. These studies suggest that the range of mechanisms by which p97 activity is controlled is vastly underexplored with significant advances possible for understanding p97 regulation by the most known adapters.

Inhibition of talin-induced integrin activation by a double-hit stapled peptide.

Integrins are ubiquitously expressed cell-adhesion proteins. Activation of integrins is triggered by talin through an inside-out signaling pathway, which can be driven by RAP1-interacting adaptor molecule (RIAM) through its interaction with talin at two distinct sites. A helical talin-binding segment (TBS) in RIAM interacts with both sites in talin, leading to integrin activation. The bispecificity inspires a "double-hit" strategy for inhibiting talin-induced integrin activation. We designed an experimental peptidomimetic inhibitor, S-TBS, derived from TBS and containing a molecular staple, which leads to stronger binding to talin and inhibition of talin:integrin interaction. The crystallographic study validates that S-TBS binds to the talin rod through the same interface as TBS. Moreover, the helical S-TBS exhibits excellent cell permeability and effectively suppresses integrin activation in cells in a talin-dependent manner. Our results shed light on a new class of integrin inhibitors and a novel approach to design multi-specific peptidomimetic inhibitors.

[Advances in molecular function of p62 protein and its role in diseases].

Sequestosome 1 (SQSTM1/p62) is a selective autophagy adaptor protein that plays an important role in the clearance of proteins to be degraded as well as in the maintenance of cellular proteostasis. p62 protein has multiple functional domains, which interact with several downstream proteins to precisely regulate multiple signaling pathways, thereby linking p62 to oxidative defense systems, inflammatory responses and nutrient sensing. Studies have shown that mutation or abnormal expression of p62 is closely related to the occurrence and development of various diseases, including neurodegenerative diseases, tumors, infectious diseases, genetic diseases and chronic diseases. This review summarizes the structural features and molecular functions of p62. Moreover, we systematically introduce its multiple functions in protein homeostasis and regulation of signaling pathways. Furthermore, the complexity and versatility of p62 in the occurrence and development of diseases are summarized, with the aim to provide a reference for understanding the function of p62 protein and facilitating related disease research.

Structural basis of ELKS/Rab6B interaction and its role in vesicle capturing enhanced by liquid-liquid phase separation.

ELKS proteins play a key role in organizing intracellular vesicle trafficking and targeting in both neurons and non-neuronal cells. While it is known that ELKS interacts with the vesicular traffic regulator, the Rab6 GTPase, the molecular basis governing ELKS-mediated trafficking of Rab6-coated vesicles, has remained unclear. In this study, we solved the Rab6B structure in complex with the Rab6-binding domain of ELKS1, revealing that a C-terminal segment of ELKS1 forms a helical hairpin to recognize Rab6B through a unique binding mode. We further showed that liquid-liquid phase separation (LLPS) of ELKS1 allows it to compete with other Rab6 effectors for binding to Rab6B and accumulate Rab6B-coated liposomes to the protein condensate formed by ELKS1. We also found that the ELKS1 condensate recruits Rab6B-coated vesicles to vesicle-releasing sites and promotes vesicle exocytosis. Together, our structural, biochemical, and cellular analyses suggest that ELKS1, via the LLPS-enhanced interaction with Rab6, captures Rab6-coated vesicles from the cargo transport machine for efficient vesicle release at exocytotic sites. These findings shed new light on the understanding of spatiotemporal regulation of vesicle trafficking through the interplay between membranous structures and membraneless condensates.

HspB8 interacts with BAG3 in a "native-like" conformation forming a complex that displays chaperone-like activity.

The HspB8-BAG3 complex plays an important role in the protein quality control acting alone or within multi-components complexes. To clarify the mechanism underlying its activity, in this work we used biochemical and biophysical approaches to study the tendency of both proteins to auto-assemble and to form the complex. Solubility and Thioflavin T assays, Fourier transform infrared spectroscopy and atomic force microscopy analyses clearly showed the tendency of HspB8 to self-assemble at high concentration and to form oligomers in a "native-like" conformation; otherwise, BAG3 aggregates poorly. Noteworthy, also HspB8 and BAG3 associate in a "native-like" conformation, forming a stable complex. Furthermore, the high difference between dissociation constant values of HspB8-HspB8 interaction with respect to the binding to BAG3 obtained by surface plasmon resonance confirms that HspB8 is an obligated partner of BAG3 in vivo. Lastly, both proteins alone or in the complex are able to bind and affect the aggregation of the Josephin domain, the structured domain that triggers the ataxin-3 fibrillation. In particular, the complex displayed higher activity than HspB8 alone. All this considered, we can assert that the two proteins form a stable assembly with chaperone-like activity that could contribute to the physiological role of the complex in vivo.

Exploration of O-GlcNAc transferase glycosylation sites reveals a target sequence compositional bias.

O-GlcNAc transferase (OGT) is an essential glycosylating enzyme that catalyzes the addition of N-acetylglucosamine to serine or threonine residues of nuclear and cytoplasmic proteins. The enzyme glycosylates a broad range of peptide sequences and the prediction of glycosylation sites has proven challenging. The lack of an experimentally verified set of polypeptide sequences that are not glycosylated by OGT has made prediction of legitimate glycosylation sites more difficult. Here, we tested a number of intrinsically disordered protein regions as substrates of OGT to establish a set of sequences that are not glycosylated by OGT. The negative data set suggests an amino acid compositional bias for OGT targets. This compositional bias was validated by modifying the amino acid composition of the protein fused in sarcoma (FUS) to enhance glycosylation. NMR experiments demonstrate that the tetratricopeptide repeat region of OGT can bind FUS and that glycosylation-promoting mutations enhance binding. These results provide evidence that the tetratricopeptide repeat region recognizes disordered segments of substrates with particular compositions to promote glycosylation, providing insight into the broad specificity of OGT.

Unique interface and dynamics of the complex of HSP90 with a specialized cochaperone AIPL1.

Photoreceptor phosphodiesterase PDE6 is central for visual signal transduction. Maturation of PDE6 depends on a specialized chaperone complex of HSP90 with aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1). Disruption of PDE6 maturation underlies a severe form of retina degeneration. Here, we report a 3.9 Å cryoelectron microscopy (cryo-EM) structure of the complex of HSP90 with AIPL1. This structure reveals a unique interaction of the FK506-binding protein (FKBP)-like domain of AIPL1 with HSP90 at its dimer interface. Unusually, the N terminus AIPL1 inserts into the HSP90 lumen in a manner that was observed previously for HSP90 clients. Deletion of the 7 N-terminal residues of AIPL1 decreased its ability to cochaperone PDE6. Multi-body refinement of the cryo-EM data indicated large swing-like movements of AIPL1-FKBP. Modeling the complex of HSP90 with AIPL1 using crosslinking constraints indicated proximity of the mobile tetratricopeptide repeat (TPR) domain with the C-terminal domain of HSP90. Our study establishes a framework for future structural studies of PDE6 maturation.

Structural Insights Uncover the Specific Phosphoinositide Recognition by the PH1 Domain of Arap3.

Arap3, a dual GTPase-activating protein (GAP) for the small GTPases Arf6 and RhoA, plays key roles in regulating a wide range of biological processes, including cancer cell invasion and metastasis. It is known that Arap3 is a PI3K effector that can bind directly to PI(3,4,5)P3, and the PI(3,4,5)P3-mediated plasma membrane recruitment is crucial for its function. However, the molecular mechanism of how the protein recognizes PI(3,4,5)P3 remains unclear. Here, using liposome pull-down and surface plasmon resonance (SPR) analysis, we found that the N-terminal first pleckstrin homology (PH) domain (Arap3-PH1) can interact with PI(3,4,5)P3 and, with lower affinity, with PI(4,5)P2. To understand how Arap3-PH1 and phosphoinositide (PIP) lipids interact, we solved the crystal structure of the Arap3-PH1 in the apo form and complex with diC4-PI(3,4,5)P3. We also characterized the interactions of Arap3-PH1 with diC4-PI(3,4,5)P3 and diC4-PI(4,5)P2 in solution by nuclear magnetic resonance (NMR) spectroscopy. Furthermore, we found overexpression of Arap3 could inhibit breast cancer cell invasion in vitro, and the PIPs-binding ability of the PH1 domain is essential for this function.

Coupled folding-upon-binding of human tumor suppressor MIG6 to lung cancer EGFR kinase domain and molecular trimming/stapling of MIG6-derived ß-hairpins to target the coupling event.

Human epidermal growth factor receptor (EGFR) is involved in strong association with malignant proliferation, which has been shown to play a central role in the development and progression of non-small cell lung cancer and other solid tumors. The tumor-suppressor protein MIG6 is a negative regulator of EGFR kinase activity by binding at the activation interface of asymmetric dimer of EGFR kinase domain to disrupt EGFR dimerization and then inactivate the kinase. The protein adopts two discrete fragments 1 and 2 to directly interact with EGFR. It is revealed that the MIG6 fragment 2 is intrinsically disordered in free unbound state, but would fold into a well-structured ß-hairpin when binding to EGFR, thus characterized by a so-called coupled folding-upon-binding process, which can be regarded as a compromise between favorable direct readout and unfavorable indirect readout. Here, a 23-mer F2P peptide was derived from MIG6 fragment 2, trimmed into a 17-mer tF2P peptide that contains the binding hotspot region of the fragment 2, and then constrained with an ordered hairpin conformation in free unbound state by disulfide stapling, finally resulting in a rationally stapled/trimmed stF2P peptide that largely minimizes the unfavorable indirect readout effect upon its binding to EGFR kinase domain, with affinity improved considerably upon the trimming and stapling/trimming. These rationally designed ß-hairpin peptides may be further exploited as potent anti-lung cancer agents to target the activation event of EGFR dimerization.

Discovery of a small molecule ligand of FRS2 that inhibits invasion and tumor growth.

PURPOSE: Aberrant activation of the fibroblast growth factor receptor (FGFR) family of receptor tyrosine kinases drives oncogenic signaling through its proximal adaptor protein FRS2. Precise disruption of this disease-causing signal transmission in metastatic cancers could stall tumor growth and progression. The purpose of this study was to identify a small molecule ligand of FRS2 to interrupt oncogenic signal transmission from activated FGFRs. METHODS: We used pharmacophore-based computational screening to identify potential small molecule ligands of the PTB domain of FRS2, which couples FRS2 to FGFRs. We confirmed PTB domain binding of molecules identified with biophysical binding assays and validated compound activity in cell-based functional assays in vitro and in an ovarian cancer model in vivo. We used thermal proteome profiling to identify potential off-targets of the lead compound. RESULTS: We describe a small molecule ligand of the PTB domain of FRS2 that prevents FRS2 activation and interrupts FGFR signaling. This PTB-domain ligand displays on-target activity in cells and stalls FGFR-dependent matrix invasion in various cancer models. The small molecule ligand is detectable in the serum of mice at the effective concentration for prolonged time and reduces growth of the ovarian cancer model in vivo. Using thermal proteome profiling, we furthermore identified potential off-targets of the lead compound that will guide further compound refinement and drug development. CONCLUSIONS: Our results illustrate a phenotype-guided drug discovery strategy that identified a novel mechanism to repress FGFR-driven invasiveness and growth in human cancers. The here identified bioactive leads targeting FGF signaling and cell dissemination provide a novel structural basis for further development as a tumor agnostic strategy to repress FGFR- and FRS2-driven tumors.

Molecular glue CELMoD compounds are regulators of cereblon conformation.

Cereblon (CRBN) is a ubiquitin ligase (E3) substrate receptor protein co-opted by CRBN E3 ligase modulatory drug (CELMoD) agents that target therapeutically relevant proteins for degradation. Prior crystallographic studies defined the drug-binding site within CRBN's thalidomide-binding domain (TBD), but the allostery of drug-induced neosubstrate binding remains unclear. We performed cryo-electron microscopy analyses of the DNA damage-binding protein 1 (DDB1)-CRBN apo complex and compared these structures with DDB1-CRBN in the presence of CELMoD compounds alone and complexed with neosubstrates. Association of CELMoD compounds to the TBD is necessary and sufficient for triggering CRBN allosteric rearrangement from an open conformation to the canonical closed conformation. The neosubstrate Ikaros only stably associates with the closed CRBN conformation, illustrating the importance of allostery for CELMoD compound efficacy and informing structure-guided design strategies to improve therapeutic efficacy.

In silico prediction, characterization, docking studies and molecular dynamics simulation of human p97 in complex with p37 cofactor.

BACKGROUND: The AAA + ATPase p97 is an essential unfoldase/segragase involved in a multitude of cellular processes. It functions as a molecular machine critical for protein homeostasis, homotypic membrane fusion events and organelle biogenesis during mitosis in which it acts in concert with cofactors p47 and p37. Cofactors assist p97 in extracting and unfolding protein substrates through ATP hydrolysis. In contrast to other p97's cofactors, p37 uniquely increases the ATPase activity of p97. Disease-causing mutations in p97, including mutations that cause neurodegenerative diseases, increase cofactor association with its N-domain, ATPase activity and improper substrate processing. Upregulation of p97 has also been observed in various cancers. This study aims towards the characterization of the protein-protein interaction between p97 and p37 at the atomic level. We defined the interacting residues in p97 and p37. The knowledge will facilitate the design of unique small molecules inhibiting this interaction with insights into cancer therapy and drug design. RESULTS: The homology model of human p37 UBX domain was built from the X-ray crystal structure of p47 C-terminus from rat (PDB code:1S3S, G) as a template and assessed by model validation analysis. According to the HDOCK, HAWKDOCK, MM-GBSA binding free energy calculations and Arpeggio, we found that there are several hydrophobic and two hydrogen-bonding interactions between p37 UBX and p97 N-D1 domain. Residues of p37 UBX predicted to be involved in the interactions with p97 N-D1 domain interface are highly conserved among UBX cofactors. CONCLUSION: This study provides a reliable structural insight into the p37-p97 complex binding sites at the atomic level though molecular docking coupled with molecular dynamics simulation. This can guide the rational design of small molecule drugs for inhibiting mutant p97 activity.

The First Class of Small Molecules Potently Disrupting the YAP-TEAD Interaction by Direct Competition.

Inhibition of the YAP-TEAD protein-protein interaction is an attractive therapeutic concept under intense investigation with the objective to treat cancers associated with a dysregulation of the Hippo pathway. However, owing to the very extended surface of interaction of the two proteins, the identification of small drug-like molecules able to efficiently prevent YAP from binding to TEAD by direct competition has been elusive so far. We disclose here the discovery of the first class of small molecules potently inhibiting the YAP-TEAD interaction by binding at one of the main interaction sites of YAP at the surface of TEAD. These inhibitors, providing a path forward to pharmacological intervention in the Hippo pathway, evolved from a weakly active virtual screening hit advanced to high potency by structure-based design.

Mechanistic Insights into the Long-range Allosteric Regulation of KRAS Via Neurofibromatosis Type 1 (NF1) Scaffold Upon SPRED1 Loading.

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

Computational studies with flavonoids and terpenoids as BRPF1 inhibitors: in silico biological activity prediction, molecular docking, molecular dynamics simulations, MM/PBSA calculations.

The BRPF1 protein is encoded by the BRPF1 gene. In addition, the BRPF1 gene is known to be upregulated in leukaemia. Recent studies have shown that it is also overexpressed in hepatocellular carcinoma (HCC) as well. Therefore, BRPF1 is a significant target for anti-cancer drug development studies, especially on HCC. 40 terpenoids and flavonoids were chosen because of their anticancer properties given in the literature. In this study, the biological activity of molecules was also investigated with in silico structure-activity relationship analysis. In addition, interactions between a series of terpenoids and flavonoids and the BRPF1 protein were investigated by molecular docking and molecular dynamics simulations. The energy change caused by the interactions of BRPF1 with different compounds was also evaluated by MM/PBSA calculations. It has been revealed that compound 5 (-9.2 kcal/mol), a kind of secoclerodane type diterpenoid, has a higher affinity both compared to other flavonoids and terpenoids, and 9F9 (-7.9 kcal/mol), a selective BRPF1 inhibitor. The study presented in this article demonstrates that compound 5, as a natural product, could form a chemical scaffold for the development of selective BRPF1 bromodomain inhibitors.